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It is shown that the dynamics of a two-dimensional crystal with a finite concentration of dislo- 
cations, as well as vacancy and interstitial defects, is governed by the hydrodynamic equations of 
a viscoelastic medium. At the longest length scales the viscoelasticity is described by the simplest 
Maxwell model, whose shear and compressional relaxation times are obtained in terms of microscopic 
quantities, including the density of free dislocations. At short length scales, bond orientational order 
effects become important and lead to wavevector dependent corrections to the relaxation times. 
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I. INTRODUCTION 



Experimental ]1|-|4| and numerical studies of the 
dynamics of driven disordered media, such as vortex ar- 
rays in type-II superconductors and charge density waves 
in metals, have clearly indicated that when the disor- 
der is strong these systems exhibit a spatially inhomo- 
geneous plastic response upon depinning, without long 
wavelength elastic restoring forces [[to). In this plas- 
tic regime dislocations proliferate due to the interplay 
of drive, disorder and interactions, and the system is 
broken up in fluid-like regions flowing around solid re- 
gions [|| . It was proposed recently |ll[] that a description 
of shear deformations in this plastic regime may be ob- 
tained by focusing on the dynamics of coarse-grained de- 
grees of freedom (the solid-like regions) and replacing the 
elastic couplings of local displacements with viscoelastic 
couplings of local velocities. It was found that the sim- 
plest, Maxwell, model of viscoelasticity yields interesting 
features, such hysteretic depinning and switching, which 
have also been observed in experiments. The local vis- 
coelastic couplings were introduced as an ad hoc way of 
mimicking the presence in the system at any given time 
of unbound dislocations, whose motion is in turn respon- 
sible for plastic slip. 

In this paper we provide some justification to the 
model studied in Ref. jllj by showing that the dynamics 
of a two-dimensional crystal with a finite concentration of 
annealed free dislocations, as well as vacancy and intersti- 
tial defects, is governed by the hydrodynamic equations 
of a viscoelastic medium (with hexatic order). Starting 
with the hydrodynamic equations for a two-dimensional 
solid with finite concentrations of dislocations, vacancies 
and interstitials obtained some time ago by Zippelius, 
Halperin and Nelson jl2| (ZHN), we show that such equa- 
tions can be recast in the form of hydrodynamic equa- 
tions for a viscoelastic medium, with a Maxwell form for 
the nonlocal transport coefficients fl3|] . 

It has of course long been recognized that the creation 
and motion of dislocations are the main mechanism by 
which a crystal undergoes plastic deformations E3,tt5[. 



A number of efforts (T^Jl^] have been made to describe 
the motion of dislocations through crystalline solids and 
relate the plastic strain rate to the dislocation dynamics. 
There has also been work specifically on the effect of dis- 
location motion on linear and nonlinear stress relaxation 



in a two-dimensional crystal |16|, where it was shown for 
a particular geometry that the crystal responds viscoelas- 
tically when free dislocations are present. Our work gen- 
eralizes this by obtaining equations that describe the vis- 
coelastic response of a two-dimensional solid to stresses 
in an arbitrary geometry. The equations incorporate vis- 
coelastic effects in the response to both shear and com- 
pressional deformations, and have precisely the structure 
of the viscoelastic equations for a simple viscous fluid 
|l8| . Furthermore, our derivation yields expressions for 
the (Maxwell) relaxation times for shear and compression 
in terms of microscopic parameters. 

Of course, in the presence of quenched disorder and 
drive, dislocations are continuously generated and healed 
in the system. Furthermore, they are not always free to 
move and relax a local stress as they may be pinned by 
disorder. For these reasons one may question the assump- 



tion behind the model introduced in Ref. |11|, namely 
that the plastic response of extended media driven over 
strong quenched disorder at zero temperature may be de- 
scribed at large scale by the same equation that govern 
the stress relaxation due to the motion of annealed dis- 
locations. At best, the effective density of "free" disloca- 
tions that can relax local stresses will be a strong function 
of the applied driving force, disorder, and possibly time. 
As shown in Ref. [pi]] , the depinning of a driven viscoelas- 
tic medium does, however, exhibit a number of features 
seen in experiments, namely memory effects and coexis- 
tence of weakly and strongly pinned degrees of freedom. 
It does, therefore provide a useful starting model for the 
description of the complex dynamics of these systems. 

Finally, it should be pointed out that crystals with an- 
nealed dislocations and simple viscous fluids are not the 
only examples of viscoelastic media. Linear viscoelastic- 
ity is of course one of the distinctive properties of com- 
plex and polymeric fluids, that can fill containers of any 
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shape and yet may shrink like rubber when stretched and 
released quickly. Molecular theories of viscoelasticity for 
complex fluids have been developed by various methods 
and there is a vast literature on the subject |19). 



II. SUMMARY OF MAIN RESULTS AND 
OUTLINE 



It is well known that a dense fluid has insufficient time 
to flow in response to a high frequency strain rate, but 
instead reacts elastically, as a solid would |^0| . This leads 
to the appearance in the fluid of propagating shear waves, 
with an associated peak at nonzero frequency in the spec- 
trum of transverse current fluctuations. Such viscoelastic 
effects are easily incorporated phenomenologically into 
the hydrodynamic description of fluids |l3[ . If a shear- 
ing force is applied to a fluid yielding a stress a xy , the 
local strain at each point can be expressed in terms of 
derivatives of the displacement field u at that point. In a 
steady-state situation the flow is purely viscous and the 
stress is proportional to the local rate of strain, 



'xy 
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du T 



dt V dy dx 



(1) 



where 77 is the shear viscosity. By contrast, if the force 
is applied suddenly, the instantaneous displacement is 
related to the stress via a typical elastic stress-strain re- 
lation, 



(du x 



8Uy 

dx 



(2) 



where (i is the instantaneous (high-frequency) shear mod- 
ulus. An interpolation between these two forms yields the 
Maxwell model of viscoelasticity, where the stress-strain 
rate relation takes the form of a differential equation, 



da xy 

dt tm ^ dy dx 



(3) 



The local flow velocity, v, of the constituent particles 
is identified as the rate of change of the displacement, 
v = dtu, and tm — r}/^ is called the Maxwell relax- 
ation time. The Maxwell model of viscoelasticity given 
by Eq. (^|) yields elastic response to high frequency shear- 
ing forces (lotm >> 1) and viscous response at low fre- 
quency (uotm << 1). 

In general, viscoelastic effects are present in the re- 
sponse to both shear and compression and can be de- 
scribed by modifying the hydrodynamic equations of a 
fluid by introducing nonlocal (in space and time) ther- 
modynamic quantities and transport coefficients. This 
can be done phenomenologically through the use of fre- 
quency sum rules as constraints, or can be formulated 



more systematically by deriving such equations via for- 
mal projection operator methods or kinetic theory tech- 
niques In its simplest realization this approach 
yields the Maxwell model of viscoelasticity. 

In this paper we show that the hydrodynamic equa- 
tions for a two-dimensional solid with a finite concentra- 
tion of free dislocations proposed by ZHN are equivalent 
to a set of hydrodynamic equations for a simple fluid 
(with hexatic order), with a Maxwell- type stress-rate of 
strain relation. For convenience, we define a symmetrized 
as 



strain rate, Vij, 



Vij = (diVj + djVi)/2 
= Vij + -o i3 v kk , 



(4) 



where we have separated the traceless part and the 
trace v kk = V • v. It is also useful to separate the stress 
tensor cr^ into its traceless part and its trace, i.e., 
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(5) 



We find that, at long wavelengths, the traceless part of 
the stress tensor is related to the strain rate according to 



~dt 



(G) 



where t 8 is the shear relaxation time, which is obtained 
in terms of microscopic parameters (see Eq. ([55]) below) . 
Equation (||) is precisely of the Maxwell form. If the 
product jj,T s is identified with the shear viscosity rj of the 
fluid, Eq. ([|) yields viscous flow in response to low fre- 
quency shearing forces. It also describes the fact that a 
liquid always flows in response to a steady shear. 

In contrast, even a liquid has a finite low frequency 
compressional modulus, that is, it resists elastically low 
frequency compressional deformations. The trace of the 
stress tensor, a k k, contains both reversible contributions 
proportional to pressure changes arising (assuming tem- 
perature fluctuations are negligible) from the fractional 
change in volume and irreversible and dissipative con- 
tributions proportional to velocity gradients. It can be 
written as 



Cfcfc 
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with 5n the deviation of the density from its equilibrium 
value n , and -Bo the low frequency bulk modulus of the 
solid (which generally depends on the density of vacan- 
cies and interstitials defects). To a good approximation 
Bq is also the compressibility of the liquid phase . We 
find that the dissipative part a kk s of the stress tensor is 



related to the local rate of strain via the equation 

Sn 



dt 



dis 
7 kk 

n 



-2(-Boo — i?o) 



where is the high frequency compressional modulus 
of the solid in the absence of defects. At high frequency 
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defects do not have time to respond to an applied force 
and therefore do not contribute to the elastic constants. 
Finally, Tb is the compressional relaxation time, which 
again is given below in terms of microscopic quantities. 
Both relaxation times, r s and Tb are inversely propor- 
tional to the density of free dislocations, u/. As a result, 
they both diverge as the melting temperature T m of the 
two-dimensional solid is approached from above as [plp2| 



, 6 cxexp(6/ (T-T m ) a ) 



(9) 



with a — 0.36963. At the longest length scales, the bulk 
viscoelasticity described by Eq. (||) is due solely to dislo- 
cation climb, while the shear viscoelasticity is due to both 
dislocation climb and glide. Dislocation climb which re- 
quires the creation of vacancies and interstitials, proceeds 
more slowly than dislocation glide which does not involve 
a change in particle number. For this reason, Tb 3> r s , in 
general. 

Of course, the 2d solid does not melt directly into a 
isotropic liquid at T m , rather it enters the hexatic phase, 
where it has only short-ranged translational order, but re- 
tains quasi-long-ranged bond orientational order (BOO). 
At the longest length scales, the presence of BOO does 
not affect the stress-rate of strain relations. When terms 
of higher order in the gradients of the hydrodynamic 
fields are included, the existence of BOO becomes appar- 
ent and modifies the stress-strain rate relations at finite 
length scales, as discussed in Section below. In particu- 
lar, this gives wavevector-dependent corrections to r s and 
Tb, corresponding to a distribution of relaxation times. 

The rest of the paper is organized as follows. In Sec- 
tion III we review the hydrodynamic equations for a two- 
dimensional solid with vacancy and interstitial defects 
[1121 . In Section IV the viscoelastic stress-rate of strain re- 



lations are derived in the longest length scale regime, and 
expressions for the long-wavelength viscoelastic memory 
times are presented. In Section |V| the effects of incorpo- 
rating higher order gradients are investigated, including 
the wavevector dependent corrections to t s . 



defect density, n A (r, t), i.e., density of interstitials minus 
density of vacancies, and the equation of motion for the 
displacement field, u(r,t), can be written as 

On 
~dt 
% 
dt 
dnA 
~dt 
du 
~dt 

where m is the mass of the constituent particles and no 
is the average number density. 

These conservation laws must be completed with con- 
stitutive equations for the two fluxes, the momentum flux 
or stress tensor, cry (r, t), and the defect current, j A (r, t). 
Assuming purely relaxational dynamics for the defects, 
j A is given by 



= --v g , 


(10a) 
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- = -v-j A , 


(10c) 
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where T is a constant kinetic coefficient and pd is the 
defect pressure, 



Pd = X 



n 



7V ■ u , 



(12) 



with x a susceptibility proportional to the sum of the 
equilibrium densities of vacancies and interstitials, and 7 
a constant associated with the coupling between changes 
in defect density and bulk distortions in the lattice. The 
constitutive relation for the stress tensor, cry, is 



cry = 2j«u/y + \w s kk S i: j +"f{Sn A /n )Sij 



(13) 



where wfj = (wy + Wji)/2 is the symmetric part of the 
strain tensor , given by 



(14) 



III. ZHN HYDRODYNAMICS OF A SOLID WITH 
DISLOCATIONS 

We begin by reviewing the hydrodynamics of the solid 
phase, following Refs. Jl2], |2j| and ]24| . We consider a 
two-dimensional hexagonal crystal with a finite density 
of vacancies and interstitials. As pointed out by Mar- 
tin |23| ] , incorporating a finite density of point defects is 
crucial to properly distinguish between local lattice de- 
formations and local changes in the number density. For 
simplicity, temperature fluctuations are neglected in all 
of the following. 

In the absence of dislocations, the equations for con- 
servation of density, n(r,t), momentum, g(r,t), and net 



The coefficients /i and A are the Lame constants. We 
note that the equilibrium bulk modulus, Bq, of the solid 
is defined in terms of the response to a static isotropic 



compression as Bq = 



J_SA 

A 8p 



, with A the area of the 



crystal. In the presence of defects it is then given by 
M + A - 7 2 X 



Bn = 



1 + (u + A + 2 7 ) X 



(15) 



On the other hand, if one were to measure the bulk modu- 
lus on time scales long compared to the elastic relaxation 
time, but short compared to the defect diffusion time, the 
effective bulk modulus, denoted by B^, would be 



Bo 



M + A 



(16) 
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corresponding to the purely elastic contribution obtained 
in the absence of defects. In a finite-frequency measure- 
ments, the value B^ > Bq is reached at sufficiently high 
frequency. Finally changes in the defect density can be 
related to changes in the overall density and bulk strains 
in the lattice, as 



Sn = Sn/ 



n w k k 



(17) 



where Sn = n — no and 5n& — ha — are the deviations 
of the density and defect density from their equilibrium 
values, respectively. 

In the presence of dislocations, it is no longer possible 
to define a continuous single-valued displacement field, 
u. One can define a single- valued, but discontinuous dis- 
placement field by introducing cuts at the location of the 
dislocations. The strain field u)y = diUj is still single 
valued and continuous, except at the location of the dis- 
locations, and it satisfies 



(18) 



where ao is the lattice spacing, tij the antisymmetric unit 
tensor, and B(r, t) is the Burgers-vector charge density. 
For a set of M free dislocations at discrete points {Ri,}, 
^=1,2, M , with Burgers vector b„, it is given by 



M 



B(r,t)=^b„J(i 



R,, 



(19) 



Dislocations are assumed to be annealed and free to equi- 
librate. The Burgers vector charge density is an extra 
hydrodynamic variable with a conservation law, 



dBj 
dt 



(20) 



where Jj is the Burgers current of the ith component 
of dislocation charge in the jth direction. For length 
scales longer than the average distance between free dis- 
locations, £d ~ l/y / riy, and time scales longer than the 
mean collision time between free dislocations, one can 
make some simplifying approximations. Specifically, one 
can ignore the effect of dynamical polarization of bound 
dislocation pairs and also processes by which neutral dis- 
locations pairs (i.e., dislocation pairs of opposite Burgers 
vectors) are created and annihilated. It was then shown 
in Ref. E3| , that the constitutive relation for the Burgers 
current can be written as 



C jkl( a kl +PdSkl) 



D$diB k . 



(21) 



The three terms on the right hand side of Eq. ( |21| ) corre- 
spond to the motion of dislocations in response to stresses 
in the lattice (the usual Peach-Koehler force), gradients 
in defect pressure and gradients in the dislocation den- 
sity, respectively. The tensors Cj fci 
by 



and D*^ are given 



Jkl 4k B T 
D„ - D, 



D 



ik _ ± ^ c [s ij s M + SjkSu) + D9 + iDc 5ji8 lk , 



(22a) 

k , 
(22b) 



where D g and D c are the diffusions constants for dislo- 
cation glide and climb, respectively. 

Th e co nserva tion laws for density and momentum, 
Eqs. ( 10a ) and ( 10b ), remain unchanged in the presence 
of dislocations. A new term must, however, be added 
to the equation for the conservation of defect density, 
reflecting the fact that a climbing dislocation acts as a 
source or sink of vacancies and interstitials. Equation 



(10c) is then replaced by 



= -Vj +n a ei k Jl 



(23) 



Since u is multi- valued, it is more convenient to consider 



an equation of motion for the strain rather than Eq. ( lOd ) 
for the displacement field. This is given by 



dw 



—dij 
n 



aae lk Jl 



(24) 



where, to linear order, g = mnov. The first two terms 
are simply obtained by differentiating the right hand side 
of Eq. (lOd). The third term arises from the presence of 
dislocations. 



The conservation laws of density, Eq. (10a), momen- 
tum density, Eq. (10b), defect density, Eq. (p3[), and 
Burgers vector charge, Eq. (]2(]), together with the con- 
stitutive equations for the currents, Eqs. (|ll|) , (|l3|) and 
(pl|), and the equation of motion for the local strains, 
Eq. (p4|) , constitute a closed set of equations to describe 
the dynamics of a 2d solid with an equilibrium concen- 
tration of free annealed dislocations and of vacancy and 
interstitial defects. 

The work by ZHN focused on the transverse and longi- 
tudinal hydrodynamic modes of these equations. These 
authors showed that the low frequency modes of a solid 
with dislocations are the same as those of the hexatic 
phase. In contrast, our goal here is to eliminate the dis- 
locations and defects degrees of freedom and obtain a 
closed set of hydrodynamic equations for the density, mo- 
mentum density and bond angle. By this procedure we 
show that the hydrodynamic equations of a solid with free 
dislocations arc equivalent, at long wavelengths, to the 
hydrodynamic equations of a viscoelastic medium with 
bond angle order. The latter consists of conservation 
laws for the density and momentum density, an equation 
of motion for the bond angle, with a constitutive equation 
for the stress tensor that takes the form of a differential 
equation. 
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IV. HYDRODYNAMICS OF A VISCOELASTIC 
MEDIUM 

Our strategy is to eliminate the defect density, tia, and 
the Burgers- vector charge density, B, from the hydrody- 
namic equations of a solid with dislocations discussed in 
the previous section. The outcome will be a closed set 
of equations for density, momentum density, bond angle 
field and stress tensor. 

It is convenient to write the strain tensor Wij as the 
sum of its symmetric and antisymmetric parts, 



Uij = (ivy + Wji)/2 + (wij - w 0i )/2 



(25) 



Only the symmetric part enters in the expression for the 
stress tensor, according to Eq. ([l3|). One can therefore 
invert Eq. ([l3]) to express wfj as a function of the stress 
tensor and of the density, as 

1 „ A + 7 



"7w = — n-,, 6 

-6, 



" 2/T y " v 4 M ( M + A + 7 ) 



7 



2(/j + A + 7) n ' 



(26) 



In order to eliminate the defects and dislocation de- 
grees of freedom, we now insert the constitutive equations 
for the defect and dislocation fluxes on the right hand 
side of Eqs. ( p9| ) and ( |30| ) and use 5n& = Sn + n w kk and 
Bj = £ikdi{w s k j + iujL)/ao to eliminate the defect den- 
sity and dislocation density in favor of density and strain 
field. Finally, the strain field is related to the stress ten- 
sor and to the bond angle field using Eqs. ( pq 1 ) and (pTj). 
This completes the transformation of the equation for 
the strain tensor into two closed coupled equations for 
the stress tensor and the bond angle field. 

For clarity, we present in this section the equations 
obtained when only terms of lowest order in the hydro- 
dynamic fields are retained. In the next section we will 
include higher order terms and show that these yield 
wavevector dependent relaxation times for shear and 
compression. At long wavelengths, the gradients of the 
defect current density can be neglected both in the equa- 
tion for the bond angle, Eq. (]2S|), and in the equation for 
the stress tensor, Eq. (|30|). Similarly, one can drop the 
contribution to the dislocation current due to diffusion 
of dislocations and approximate 



J k 



- C LMnl +PdSnl) 



(31) 



where Eq. (|17]) was used to eliminate the defect density 
in favor of the density and the strain. The antisymmetric 
part, wfj, has only one independent component which is 
simply the bond angle field, 9(r,t) p5| , 

(27) 

The equation of motion for the strain, Eq. J24]), can then 
be split into an equation of motion for the symmetrized 
strain, given by 



at 



k I > 



(28) 



and an equation for the bond angle field, 



(29) 

We note that the trace of Eq. (|2^) simply gives the con- 
tinuity equation for the density. 

We can now use Eq. ( p6|) to rewrite the equation of 
motion for the symmetric part of the strain tensor as a 
differential equation for the time evolution of the stress 
tensor, with the result 



d<Ti 



+a fJ,(e ik Ji + ejkJlj + a (X + 7)%e w j\ 
n V J J n 
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X~ 1 (m + A)- 7 2 Sn 
fi + A + 7 n 



(32) 



Carrying out some tedious algebra, we obtain the desired 
set of hydrodynamic equations, consisting of exact con- 
servation laws for density and momentum, the equation 
for the bond angle field [Eq] , 



? = iz.(Vxv) + ^V 2 0, 

at 2 v ' 2 

and the constitutive equation for the stress tensor, 

daij 1 _ 1 
dt t s 2r b 

= 2^Vij + BooSijVkk -Sij — , 

n no 



(33) 



(34) 



where we recall that = /i+A. The shear and compres- 
sional relaxation times (or "Maxwell" times) are given by 



1 



D n + D r 



n 



r 



Hal 



_L_ = D c B eff x 1 o 2 o 

n ~ Uf k B T~ 



Bo 



(35) 
(36) 



where B e f f = /i + A — X7 2 is the "effective" bulk mod- 
ulus that one would obtain by integrating out the defect 
density fluctuations from the free energy. 
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(a) The dislocation begins on the left side of the 
crystal, i.e. dx=0. The initial relative displacement 
of the upper and lower halves of the crystal is zero. 




(b) The dislocation has glided two lattice spacings to 
the center of the crystal, i.e. dx=Lx/2. The relative 
displacement of the upper and lower halves of the crystal. 



(c) The dislocation has now glided all the way to the 
right of the crystal, i.e. dx=Lx. The relative displacements 
of the upper and lower halves is now one lattice spacing. 

FIG. 1. An illustration of how dislocation glide leads to 
relative displacement of two half crystals. 

The shear relaxation time has a simple physical in- 
terpretation. Considering, for simplicity, only disloca- 
tion glide, the motion of a dislocation of Burgers vec- 
tor b = x a distance dx along its slip line contributes 
an amount jf-dx to the displacement of the half crys- 
tal, with L x the length of the crystal along the slip 
direction. This is illustrated for a square lattice in 
Fig. (|l]). The corresponding change in shear strain is 
~ (aodx/ L x )/ L y , with L y the width of the crystal in 
the direction normal to the slip line. The change in 
the strain 7 in a time dt due to M dislocations is then 



M dx _ 
LzL v ^ dt 



njCLQVg, with v g 



di 



the dislocation 



glide velocity. This is the well-known Orowan law |T^] 
relating the plastic strain rate to the dislocation velocity. 
Dislocation motion is driven by local stresses and we can 



write v g — H g F s , with fx g — Dg/^ksT) the dislocation 
mobility and F s ss cr s ao ~ /xoo the force on a dislocation 
due to the shear stress, a s . If we define a shear relaxation 
time by -3; ~ and assume 7 ~ 0(1), we immediately 

obtain l/r s ~ nf-j^^fiaQ. When both glide and climb 
can occur the mobilities are simply additive. 

An alternative interpretation of Eq. ([35]) is as fol- 
lows. The rate at which dislocations diffuse thermally 
(via glide and climb) a distance of the order of the 
typical separation ~ n * between free dislocations is 
l/rth = (D g + D c )nf . This is just the mean time for col- 
lisions between free dislocations. The ratio of the rate of 
diffusion of dislocations under the influence of shear elas- 
tic forces to the rate of thermal diffusion is of the order 
of the ratio of a typical shear elastic energy, /idg to ther- 
mal energy, fc^T. With this energetic argument we esti- 
mate l/r s ~ 1/rjh (/xflg/fesT) which is the result given in 
Eq. (|3g) . Interactions between dislocations are negligible 
only at time scales longer than the mean collision time, 
Tth As a result our model applies only if r s > Tth- Since 
lia^/kB ~ 0(T m ), the condition r s > Tth simply corre- 
sponds to being above the melting temperature, T m . 

It is not surprising that the bulk memory time is in- 
dependent of Dg, since pure compressional strains can 
only relax through dislocation climb. Also, since dislo- 
cation climb involves a change in particle number, D c is, 
in general, much smaller than D g , so that t s << t&. In 
addition, above the melting temperature the response to 
shear distortions is purely viscoelastic while the response 
to bulk distortions is effectively still elastic at large scales, 
with some dissipative effects (due to the climbing of dislo- 
cations) which, due to the smallness of D c , are expected 
to be negligible. 

In the limit of nj — > 0, i.e., as T m is approached from 
above, the shear and compressional relaxation times di- 
verge, T s b — > 00 according to Eq. (^), and we obtain 

ddi 



dt 
dt 



2/AV 



2 (-Boo — B Q )vkk 



t Bqp - B dSn 
dt 



n 



(37a) 
(37b) 



In t he ab sence of dislocations Vij = dtwfj. Integrating 



Eq. (37a) then yields the shear stress-strain relation for 
an elastic solid, 



2/iK 



U s 



(38) 



Similarly, to lowest order in the gradients, we can neglect 
the rate of change of defect density. Integrating Eq. ( |37b| ) 
and using Eq. (p^), we then obtain a^ k s = (Bm — B )w s kk 
and 



CTfcfc = 2BOOW 



(39) 



Far above T m , free dislocations proliferate and both 
T,, and Tb become very small. In this limit Eq.(|3^) then 
reduces to 
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Sti 

aij = 2[ir s Vij + 6 i j(B 00 - i? )T & %V • v - SijB — . (40) 

no 

Identifying the combinations fir s with the shear viscosity, 
77s, and the combination {B^ — Bo)ti, with the bulk vis- 
cosity, 77f,, Eq. ( |4C| ) is just the stress-rate of strain relation 
of a viscous liquid. 

When only terms of lowest order in the gradients of the 
hydrodynamic fields are retained, it is not apparent from 
the viscoelastic stress-rate of strain relation that a crys- 
tal with annealed dislocations is a hexatic, rather than a 
isotropic viscous liquid. The presence of quasi-long-range 
bond orientational order in the hexatic manifests itself in 
a resistance to distortions that produce a non-zero gra- 
dient in the bond-angle field, 0. Contributions to the 
hydrodynamic equations from gradients of correspond 
to spatial gradients of Wij , which have been neglected in 
this section. Such terms must be retained to see a sig- 
nature of bond orientational order in the stress-rate of 
strain relations. This is done in the next section. 



V. INCLUSION OF HIGHER ORDER 
GRADIENTS 

In this section we incorporate all diffusive terms that 
were neglected in the previous section. We do so, how- 
ever, only for a simplified model where dislocation climb 
is neglected. This corresponds to freezing out fluctua- 
tions in the defect density. It is an excellent approxi- 
mation in most cases as the diffusion constant D c for 
dislocation climb is generally very small due to the large 
energy cost of creating vacancies and interstitials. In the 
absence of fluctuations in the defect density the relation 
between symmetrized strain and stress simplifies to 



A 



I 

2/^ j ~~ 4/x(/x- 



A) 



(41) 



The exact conservation laws for density and momentum 
are of course unchanged and are still given by Eqs. (10a) 
and (10b), respectively. The equation of motion for the 
bond angle field now contains terms due to dislocation 
diffusion, which couple to gradients of stress and it is 
given by 



g = h ■ (V X v) + ^0 ^ufkdM ■ ( 42 ) 

^From Eq. (|42|), one sees that the presence of gliding 
dislocations leads to a dissipative term oc V 2 # and a cou- 
pling tO 

In the absence of fluctuations in the defect density, the 
stress tensor's trace is simply 



Okk 



Sn 

-2Boo — 
no 



(43) 



We therefore focus on its traceless part, cr# Inclusion of 
the extra diffusive terms in the constitutive relation for 
aij will yield terms involving the gradient of the Burg- 
ers vector charge density. These terms are easily recast 
as gradients of bond angle and stress, resulting in the 
following constitutive relation for cry 



dt 



a ij D g 2 - 
— V <7j 

2fiv. 



V D g 
2 

di dj - 



{tikdjdk + e jk did k )9 



Sn 



V D 9 

4 V"'" J 2 / /;„ 



(44) 



where Eq. ( |43] ) relating akk and Sn/no has been used. It 
is instructive to rewrite Eq. (|44|) in Fourier space. The 
equation for cry (k, t) is then given by 



M (k i k j -8 ij /2)ek' 2 -, (45) 

V / Tin 



t s (0) v " - / /),, 
with ki = ki/k and where we have defined 



~VD 9 



(46) 



One is naturally led to introduce a finite wavevector shear 
relaxation time, T s (k), given by 



T.(k) 



1 + £ 2 k 2 ' 



(47) 



When the constitutive relation is written in this form, it 
is apparent that £ is the length scale beyond which the 
small gradient approximation is valid. At low frequency 
Eq. (I45|) reduces to 



2ij s _ /j,D g T s (0) 



1 + £ 2 fc 2 



2(l + £ 2 fc 2 ) 



(tikkjkk + tjkkikk) I 



/' 



1+ek 2 



Sn 
no 



(48) 



By comparing Eq. ( fig ) to the corresponding relation for 
a hexatic given in ZHN, 



~ hex 



2r\ s v v 



K, 



(49) 



it is natural to identify the combination /iD g T s (0) with 
Ka, the elastic constant associated with long wavelength 
distortions of the bond angle field. The length scale £ 
can then be written as 



(50) 
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and it represents the length scale below which the energy 
associated with bond angle distortions becomes compa- 
rable to that associated with spatially uniform shear dis- 
tortions. 

The effective shear relaxation time, r s (fc), is larger at 
short length scales, as the additional stiffness due to bond 
orientational order makes the system more solid-like. Fi- 
nally, it should be noted that the term in Eq. ( [45| ) oc 
Sn/n only contributes to the equation for the longitu- 
dinal (oc V • g) component of the momentum, but not 
to that for the transverse part (oc V x g). This indi- 
cates that, even if dislocation climb is forbidden, at short 
length scales the longitudinal component of momentum 
is affected by dislocation glide, which yields higher order 
pressure-like gradients in the equations of motion. 

Clearly a regime where the wavevector dependence of 
the shear relaxation time is important only exists pro- 
vided £ > with £d ~ the average distance 
between free dislocations. At length smaller than £d in- 
teractions among dislocations can no longer be neglected 
and our model becomes invalid. It is easy to see that 



(51) 



with Tm the melting temperature of the two-dimensional 
lattice. Since our work always describes the region 
T > Tm, where free dislocations are present in the crys- 
tal, the wavevector dependence of the relaxation time will 
generally be important. 

It is a pleasure to acknowledge discussions with Alan 
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